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POISSON LIMITS OF SUMS OF POINT PROCESSES AND A 
PARTICLE-SURVIVOR MODEL 

By Matthew O. Jones and Richard F. Serfozo 

Austin Peay State University and Georgia Institute of Technology 

We present sufficient conditions for sums of dependent point pro- 
cesses to converge in distribution to a Poisson process. This extends 
the classical result of Grigelionis [Theory Probab. Appl. 8 (1963) 172- 
182] for sums of uniformly null point processes that have Poisson 
limits. Included is an application in which a particle-survivor point 
process converges to a Poisson process. This result sheds light on 
the "surprising" Poisson limit of the species competition process of 
Durrett and Limic [Stochastic Process. Appl. 102 (2002) 301-309]. 

1. Introduction. Tilman [7] studied a model for species coexistence where 
species with higher death rates are superior competitors (as can be the case 
with parasites that kill their hosts quickly). He described the dynamics of 
such systems with differential equations and provided solutions describing 
the relative abundances of the species in equilibrium. Interestingly, he was 
able to show that inferior competitors with suitable traits can endure, re- 
gardless of the number of species in the system. 

To further examine the number of species that can coexist, May and 
Nowak [5] considered a model where new species continually arrive to a sys- 
tem, possibly due to mutations or arrivals of new species from outside the 
system. They assumed sufficient time would elapse between arrivals for the 
system to come to equilibrium. Among other things, the authors were in- 
terested in the limiting distribution of the numbers of different species with 
death rates in subintervals of (0, 1). In order to study this, they introduced 
a toy model where an arriving species with death rate v immediately annihi- 
lates any species with a smaller death rate with probability a, independently 
of everything. 
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Durrett and Limic [1] gave a formal treatment of this model. They studied 
a system in which particles arrive to the interval (0, 1) according to a Pois- 
son process, and the location of each particle in the interval has a uniform 
distribution, independent of everything else. The location of a particle rep- 
resents the particle's rank, which affects how long it will survive. Whenever 
a new particle arrives, each particle currently in the interval with a rank 
below that of the arrival is independently deleted with a fixed probability. 
The evolution over time of the point process of particles in the interval is 
represented by a Markov jump process on a space of counting measures. The 
authors proved this Markov process converges to a Poisson process on (0, 1); 
its intensity is increasing and approaches infinity at 1. 

In trying to determine why Poisson processes arise as limits in such appli- 
cations, we studied several types of sums of point processes with dependen- 
cies that are abstractions of those in the particle process. This led us to a 
general Poisson convergence theorem for sums of dependent point processes 
(Theorem 1). This result is similar in structure to Grigelionis' [2] theorem 
that gives necessary and sufficient conditions for independent uniformly null 
processes to converge to a Poisson process. 

The rest of this study shows that Theorem 1 is a natural framework for 
establishing Poisson limits of particle-survivor processes. We apply it to a 
system in which particles enter K at arbitrary random times, and their loca- 
tions in 3? are independent and identically distributed with a general, con- 
tinuous distribution. Each new arrival may delete any particle located below 
it with a location-dependent probability a(x) (see Section 3). We show that 
the point process of particles converges in distribution to a Poisson process. 
Our results indicate that, by appropriately choosing the function a(x), one 
can allow an arbitrarily large number of particles to exist in any subinterval 
of (0,1), on the average. This may appear contrary to the "limiting simi- 
larity" of Tilman's model, but is actually a simple consequence of the toy 
models of [1, 5] since each particle has a finite lifetime w.p.l. 

The contents of the remaining sections are as follows. Section 2 contains 
our main result for sums of point processes and several corollaries for special 
cases. A preliminary lemma on the convergence in distribution of products 
of random variables is of interest by itself. Section 3 describes the particle- 
survivor process described above and its Poisson limit. Sojourn times of 
particles in this system and properties of departures are the topics of Sec- 
tion 4. We end in Section 5 by giving necessary and sufficient conditions for 
the stationary distribution of a Markovian particle-survivor process to be 
that of a Poisson process. 

2. Sums of dependent point processes. Before stating our main result, 
we review some terminology on point process as in [3, 4]. Let (fi,.? 7 , P) be 
a probability space and E be a Polish space endowed with its Borel cr-field 
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£. A point process £ on E is a mapping from £1 to the set of integer- valued 
measures of the form J2k $x k (•) that are finite on compact sets, where S x (-) is 
the Dirac measure with a unit atom at x. This space of measures is endowed 
with the Borel u-field generated by the topology of vague convergence. Then 
£(-B) =J2k5x k (B) represents the number of points in B G £, where Xj- are 
the point locations. The mean measure of £ is fJ-(B) = E[£(B)], B G £. 

The point process £ is a Poisson process with mean measure fi if £(B±), . . . , 
£(B m ) are independent for disjoint B\, . . . , B m G £, and £(B) is a Poisson 
random variable with mean [J,(B), where fi is finite on compact sets in 6. 
Here £ denotes the class of relatively compact sets in £. We also let £^ 
denote the sets in £ whose boundary has /^-measure 0. 

Let C~k denote the space of nonnegative functions on E with compact 
support. The integral of a functions / G with respect to a point process 
£ is the summation 



Similarly, fj,f = f E f(x)d[i(x) for a measure //. 

A sequence of point processes £ n converges in distribution to the point 
process £, denoted by £«—►£, if the probability measure of £ n converges 
weakly to the probability measure of £. Each of the following statements 



is equivalent to £ n 4 £: (a) £„/ S £/, f G C+. (b) E[e~^] -» ^[e"^], / G 
C+. (c) (^(J3i),...,en(S fc ))-i(^(J3i),...,e(fl fc )), B 1 ,...,B k E£„. When 



the limit £ is Poisson with mean measure fi, statement (b) for the Laplace 
functional convergence is equivalent to 



We are now ready for the main result. Let £ n i,£n2> ■ • • > for n > 1, be a 
finite or infinite doubly indexed array of point processes on a Polish space 
E with Borel u-field £, so that £ n fc(-B) is the number of points in B G £ . 
Assume the processes are defined on a single probability space P), 
and there are cr-fields J- n , n > 1, in J- such that £ni,£n2> ■ ■ • are conditionally 
independent given T n . Assume that either 




k 



E[e 5n/ ] ^exp 




(1) 



J2 p tink(B)>l\T n } 2 ^0 asrwoc, Pg£ 



k 



or that 



(2) 



supP{£ nfc (£)>l|j;„}4o 



as n 



oo, Be£. 



k 
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Consider the sequence 

(3) N n = Y / Znk, n>l, 

k 

and, for the infinite sequence case, assume that N n is finite on compact sets, 
so that it is a point process on E. Let N denote a Poisson process on E 
with mean measure [i. The following result gives conditions under which 

Theorem 1. For the processes defined above, 

(4) P{N n e-\F n }^P{N E-} asrwoo, 
if and only if the following conditions hold as n — > oo : 

(5) Yl P tink(B)>2\F n } 4o, B £ £, 

k 

(6) J2 P tink(B)>l\F n } Be^ 

Moreover, these conditions imply N n — > iV, as n —> oo. 

This result is a conditional-distribution variation of the classical theorem 
of Grigelionis [2] (which is Theorem 16.18 in [4]). In particular, he proved 
Theorem 1 for independent £ n i ; £n2 ; ■ • ■ that satisfy the uniformly null prop- 
erty (2), and all the conditional distributions are ordinary ones. In this case, 

his result states that (5) and (6) are necessary and sufficient for N n N. 

In addition to the standard uniformly null property (2), we included the 
alternative property (1), which we use later. 

Our proof below of Theorem 1 uses ideas and techniques associated with 
the convergence of point processes (e.g., see Lemma 5.8 and the proofs of 
Theorems 5.7 and 16.18 in [4]). The following preliminary result plays a key 
role. 



Lemma 2. Let Y n \,Y n 2, . . . , for n > 1, be a (possibly finite) doubly in- 
dexed array of random variables in (0, c], where c < 1, and let Y be a non- 
negative random variable. Assume 

(7) Y^ k -i or supY n k —> as n — > oo. 

k k 

Then as n — > oo, 

Y[(l - Y nk ) 4. e~ Y if and only if ^ Y nk 4- Y. 
k k 
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Proof. The assertion is equivalent to 

(8) - 5>g(l - y nfc ) 4 Y if and only if ^y nfc Ay. 

k k 

To prove this equivalence, it suffices to show D n — > 0, where 

(oo 
E (XrTk~ 2 M 
m=2 

Using Y nk < c and a n = sup fe Y nk , we have 

(9) D n <(l- C y l ^2Y^ k , 



(10) D n < Y nk < E M 1 - *W 

c fc 1 c k 



Then L> n A follows from (9) under the supposition J2k^nk ~* 0- Also, 
under the alternative supposition a n — > 0, if one of the limit statements in 
(8) is true, then D n — > by (10), which proves the other limit statement in 
(8). □ 

Proof of Theorem 1. We begin by determining a convenient expres- 
sion for the statement P{N n € -{Fn} P{N £ •}. This convergence is equiv- 
alent to the random Laplace functional convergence (e.g., see Theorems 5.3 
and 16.16 of [4]) 

(11) E[e- N »f\T n ]±E[e- N f], feC+. 

Since £ni,£n2> • • • are conditionally independent given T n% 

E\e- N ^\T n \ = \{E[e-^f\T n ] = - Y nk ), 
k k 

where 

Y nk = E[l-e-^ f \T n }. 

Also, the Laplace functional of the Poisson process iV has the well-known 
form E[e~ N f] =e~^ h , where h(x) := 1 — e~^ x \ Then (11) is equivalent to 

(12) \{{l-Y nk )Ae-^\ feC+. 

k 

Keep in mind that Y nk is a function of /. 
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We will complete the proof by applying Lemma 2 to establish that (5) 
and (6) are necessary and sufficient for (12). Clearly, Y n i,Y n 2, ■■ ■ are in 
(0, 1 — e~ b ], where b = max xe E f(x). Next, note that 

(13) 1 - e-*»* f < l{Znk{S f ) >1), 

where Sf is the support of /. Then by assumption (1), 

J2Ynk<J2 P {U(S f )>l\T n } 2 S0. 

k k 
Also, by (13) and assumption (2), 

supy nfc < sup P{£ nk (Sf) > l\Fn} 0. 

k k 

In light of the last two limit statements, Lemma 2 says that (12) is equivalent 
to 



(14) ]T(1 - E[e-e*»f\F n ]) =Y,Ynk ^Vh, f € C+ 

k k 

Therefore, it remains to show that (5) and (6) are necessary and sufficient 
for (14). We show sufficiency first. 
We can write 

Y,Y n k = £ E[(l - e~^)l(U(Sf) = l)\? n ] 

, . k k 

(15) 

+ - e-^ f )l(U(S f ) > 2)\F n }. 

k 

The last sum is bounded by J2k B{Cnk(Sf) > 2)|jF n } which converges in 
distribution to by assumption (5). The first sum on the right-hand side in 

(15) equals ~£,k E [lnkh\^n] =Vn,h, where 

Lk(B) = U(B n s f )i(u(Sf) = i), 

T) n (B) = J2E[Lk(B)\Fn] = £ P tink(B fl S f ) = l\F n }. 

k k 

Now, assumptions (5) and (6) imply i] n A /x, and so rj n h — » fih. Applying 
the preceding observations to (15) proves that (5) and (6) are sufficient for 
(14). 

Now suppose (14) is true and recall that (14), (12) and (11) are equiva- 
lent. Then from (11) and switching from Laplace transforms to generating 
functions, it follows that, for any fixed B € £„, 

E[s N «W\f n ]SE[s N( - B \ se[0,l]. 
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This in turn implies that (12) and, hence, (14) also hold with e~^ nk $ replaced 
by s?«fc(- B ). In particular, the version of (14) says for s G [0, 1] that 

(16) H n {s) =Y^~ E(s^ B ^\T n )] 4 (1 - s)»{B). 

k 

Here we are using the fact that E[e~ aXn ] —> E[e~ aX ], a G [0, oo), implies 
X n — >■ X, which implies E[s Xn ] — >■ E[s x ], s G [0, 1]. Then (6) follows since 

(17) Y, P{U(B) > l\T n ) = fln(0) ^ fi(B). 

k 

Next, note that 



l-Y J S m P{i nk {B)=m\T n } 



m=0 



= (1 - s)H n (0) + E E ( s " s m )Ptink{B) = m\F n }. 

k m=2 

Then using this expression along with (16) and (17), 

(s - s 2 ) E P{Cnk{B) > 2|jF n } < H n {s) - (1 - s)H n (0) ^ 0. 

k 

These observations prove that (14) implies (5) and (6). 

This completes the proof of the first assertion that (5) and (6) are neces- 
sary and sufficient for (4) . Finally, (4) and the bounded convergence theorem 

yield N n N, which proves the second assertion. □ 

We will now consider the point processes defined above, where £ n fc are 
single-atom measures. Suppose 

(18) N n (B) = J2U nk Sx nk (B), B £ £, 

k 

where X n i,X n 2, ■ ■ ■ are E- valued random variables, and U n \,U n 2, ■ ■ ■ are ran- 
dom variables that take on values or 1. Notice that (18) is the same as (3) 
with £ n k(B) = U n kSx nk (B). Assume that, for each n > 1, there is a u-field 
T n in T that contains a(X n i,X n 2, ■ ■ ■), and U n i,U n 2,... are conditionally 
independent given T n . Define 

r nk = P{Unk = M^n}, k>l. 
Corollary 3. Suppose the N n defined in (18) are such that 



(!9) J2 r nk S x nk (-)^0 and E r nk5x nk (■) n, 



as n — > oo, 
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for some measure fi that is finite on sets in £. Then N n — > N, as n—> oo, 
where N is a Poisson process with mean measure fj,. 

Proof. Note that P{U nk S Xnk (5) = l\T n } = r nk 5 Xnk (B). Then the con- 
ditions in (19) imply (1) and (6), and (5) is trivially satisfied. Then N n N 
follows by Theorem 1. □ 

Example 4. Random variable convergence. As a special case of (18), 
suppose N n = J2kU n k, where U n \ , U n 2 , ■ ■ ■ are random variables that take 
on values or 1. Assume that, for each n > 1, there is a cr-field J- n in T 
such that U n i,U n 2, ■ ■ ■ are conditionally independent given T n . In addition, 
assume that as n — > oo , 

J2P{Unk = l\Fn} 2 ^0 and J2 P { U nk = MFn}^V, 

k k 

for some [i > 0. Then by Corollary 3, the random variables N n converge in 
distribution to a Poisson random variable with mean \i. 

The next result is a special case of Corollary 3 that we use for the particle 
system model in the next section. Suppose 

n 

(20) N n {B) = Y J U nk 5 Xk {B), Be£, 

k=l 

where X±,X2, ■■■ are independent identically distributed random elements 
in E with distribution F, and U n ±, . . . , U nn are random variables that take on 
values or 1 and are conditionally independent given J- n = o~(Xi, . . . ,X n ). 

Corollary 5. Suppose the point processes N n given by (20) satisfy 

(21) {r nk ,X k ) =(r n i,X!), 1 < k < n, 

(22) lim n [ E[r 2 nl \X x =x\dF{x) =0, Be£. 
n ^°° Jb 

In addition, assume there is a function r:E^ [0, oo) such that 

(23) lim / E[\nr nl -r(x)\\X 1 = x]dF(x) = 0, B G £, 

n -^°°JB 



and the measure fJ,(B) = f B r{x)dF(x), B € £ , is finite on compact sets 
Ti 



Then N n —*N as oo, where N is a Poisson process with mean measure 



POISSON LIMITS OF SUMS OF POINT PROCESSES 



9 



Proof. The assertion will follow by Corollary 3 upon verifying (19). 
For any / € C^, with support Sf, it follows by (21) and (22) that 



E 



Y, r nkf(Xk 



Lk=l 



<maxf(x) [ nE[rl 1 \X 1 =x]dF(x) ^0. 

Js, 



This implies J2k=i r nkf(Xk) ~* 0, which proves the first condition in (19). 
To prove the second condition in (19), it suffices to show the Li-convergence 



(24) 



E 



r nkf{X k ) ~ flf 



k=l 



0. 



fee 



K • 



By the classical L\ law of large numbers, the independence of X\,X2,. 
and E\r{X\)f {X\)\ = fxf, we have 



E 



n~ l Y,r(X k )f(X k )-[if 
k=l 



0. 



Also, by (21) and (23), 



E Y. r nkf{X k ) -n" 1 ^r(X fe )/(X fc ) 
Jc=l k=l 

< max/(a;)£'[|nr ri i — r(Xi)|] — > 0. 

X 

The last two limit statements imply (24). □ 



3. Particle system with attribute-dependent survival. Suppose particles 
(representing microbes, molecules, customers, etc.) arrive to a system at 
finite times < T\ < T2 < ■ ■ •. The number of particles that arrive in the 
time interval [0,t] is A(t) = Y^=i ^ t)i this is finite when T n — > 00 
w.p.l, which we assume for convenience. No other conditions are imposed 
on the dependency among the T n . Without loss of generality (as explained 
in Proposition 9), we assume the system is empty at time 0. 

Each particle that arrives is labeled by a real-valued attribute (or rank) 
x. We interpret x as the "location" of the particle in the attribute space, 
and refer to a particle with attribute x as an "x-particle." We let X n denote 
the attribute of the particle that arrives at time T n , and assume X n are 
independent, continuous random variables, independent of the arrival times, 
with a common distribution function F. Then the space of the particle 
attributes is the open interval 



E={xGK:0<F(x) < 1}. 
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The particles in the system (i.e., in the space E) are subject to deletion 
by future arrivals as follows. Whenever a new particle arrives, it considers 
for deletion all particles currently in the system that have strictly lower 
attributes. Specifically, an arriving y-particle considers all x-particles with 
x < y for deletion, and then it deletes each such x-particle with probability 
a(x), independently of everything else. The probability that an a>particle 
survives £ deletion attempts is therefore (1 — a(x)) e . Prom the view of an x- 
particle, the probability it is deleted by a y-particle is l(x < y)a(x). Also, the 
probability that an x-particle is considered for deletion by the next arrival 
is F(x) = 1 — F(x), and the probability it is actually deleted by that arrival 
is a(x)F(x). 

We denote the lifetime of the -X^-particle by L n . That is, L n = I means the 
X^-particle dies (exits the system) at the £th deletion attempt. Under the 
preceding assumptions, the pairs (X n , L n ), n > 1 are independent, identically 
distributed, independent of the arrival times and 

(25) P{L n > t\X n } = (1 - a(X n )f, £>1. 

The only other assumption we make is that the function l/a(x) is bounded 
on compact sets in E. This implies the measure \i defined below in (27) is 
finite on compact sets. 

We will consider the point process Nt on E, where Nt(B) denotes the 
number of particles with attributes in B that are still in the system at time 
t. That is, 

A(t) 

(26) N t (B) = J2HL k >Q tk )5 Xk (B), Be£, 

k=l 

where Qtk = J2f= k +i ^-i-^-j > Xk) 1S the quantity of attempted deletions of 
the Xfc-particle by those X,-particles that are above and arrive in the 
interval [0, t] after X^. Recall that the dynamics are such that the A^-particle 
is alive at time t if and only if its lifetime exceeds Qtk ■ 

The time-dependent evolution of the particle process {Nt : t > 0} is highly 
dependent on the nature of the arrival process A(t), which may have any 
probabilistic structure. However, when the arrival process is Poisson, then 
the particle process is a pure-jump Markov process. Specifically, the particle 
process changes state only at the Poisson arrival times T n , and the embedded 
sequence of states Nt u at the arrival times is a Markov chain. That is, Nt 
is a Markov chain subordinated to the Poisson process A(t). 

In this Markovian setting with the additional assumptions that the dele- 
tion probability a(x) = a is independent of the location x and F is a uniform 
distribution on E = (0, 1), Durrett and Limic [1] showed that Nt converges 
in distribution to a Poisson process with mean measure (J>(B) = j B a ^_ x ^ dx. 
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The approach they used to prove this is not applicable in the non-Markovian 
setting with location-dependent deletions. 

We are now ready to apply the results above for sums of point processes to 
show the particle process N t converges in distribution to a Poisson process. 

Theorem 6. Under the preceding conditions, Nt—>N as t — > oo, where 
N is a Poisson process on E with mean measure 

(27) fj.(B)= I 1 dF(x), Be£. 

jb a(x)r (x) 

Intricate dependencies in the particle process cloud the fact that the Pois- 
son limit of Nt follows by using the framework in Corollary 5. Furthermore, 
the main property producing the Poisson limit is the claim (29) below about 
ranks of order statistics (but in terms of conditional distributions). Further 
insights underlying Poisson limits are in Section 5, where we discuss station- 
ary distributions. 

Proof of Theorem 6. We begin with the key observation that, con- 
ditioned on A(t) = n, the process Nt depends on time only through the 
"order" in which the n particles arrive. Because of this observation, we can 
formulate the distribution of Nt as follows. 

Fix an n > 1, and let (X±, L±), . . . , (X n ,L n ) be i.i.d. versions of the ran- 
dom vectors above, but here the subscripts are not the order of particle 
arrivals. On the same probability space and independent of the {X k ,L k ), let 
7r n i, . . . , 7r nn be the random permutation of 1, ... ,n, where each permutation 
is equally likely with probability 1/ro!. We define the point process £ n on E 
by 

n 

(28) Zn = J2 U nk$X k , 



k=l 



where U nk = l(L k > Qnk) and 



Qnk = > X k, Knj > 7T„fc), 1 < k < n. 

3=1 

We interpret X±, . . . ,X n as the attributes of the first n particles to arrive, 
with 7r n i, . . . , TTnn representing the order in which they arrive to the system; 
for example, ix n k = 3 means that the X^-particle is the third one to arrive. 
Also, Qnk is the quantity of attempted deletions of the X^-particle. 

From the definition of the process Nt, it is clear that £ n represents Nt 
conditioned on A(t) =n. In other words, 

P{N t e-} = E[P m (.)], 
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where P n (') — P{£,n £= •}, n > 0, is the distribution of £ n , and £o is the zero 
measure. Since P n does not depend on t and ^4(t) — > oo w.p.l, the assertion 

Nt — > N will follow upon showing that -P{£ n € •} converges weakly to P{N € 

•}, or equivalently, that £ n -^ N. 

We will prove this convergence by applying Corollary 5 to the process £ n 
in (28). We begin by proving the following assertion. 

Claim. The Q n i, ■ ■ ■ , Qnn are conditionally independent given T n = o~(Xi, 
...,X n ), and 

(29) P{Qnk = m\J 7 n } = l/v n k, < m < u nk - 1, 1 < k < n, 
where v nk = YJj=i 1 ( x j >X k ). 

Proof. It suffices to prove the claim for {X\ > • • • > X n }. In this case, 
v n k = k and 

fc-i 

Qnk = ^2 1 ( 7r «i > n nk), 

which is the number of the particles X\, . . . , X k _\ that arrive after X k . Since 
7r n i, . . . , ir nn is the equally likely permutation of 1, ... ,n, we have 

P{Qnk = mk\F n } = V^nJfc) < m k < v nk - 1, 

(30) 

( n \ n 

P\ fl {Qnk = m k }\T n = 1/n! = J] V^n*- 

U=i J fc=i 

The first equality of the second line above follows because, conditioned on 
J- n , each of the n\ equally likely permutations of arrival orderings corre- 
sponds to exactly one distinct assignment of Q nk values. This proves the 
claim. □ 

Next, we derive a simple expression for r nk = P{U nk = l^n}- Using (25) 
and a(x) = 1 — a(x), we have 

r nk = E[P{U nk = l\Q nk ,Fn}\fn} = E[a{X nk ) Q ^\F n }. 

Then from (30), 

(31) r nk = 2^ a K X nk) = (v \ • 

v nk m=Q V nk a{A nk ) 

We are now ready to verify the assumptions in Corollary 5. By assump- 
tion, X nk = X k , 1 < k < n, are independent and identically distributed with 
distribution F, and U n ±, . . . , U nn are conditionally independent given T n by 
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the claim we proved above. Also, (r n k,Xk) = (r n i,X\), 1 < k < n, follows 
from (31) and the fact that v n k = Y^=i > ^k) is a function of X k and 
the set {X 1 ,...,X n }\{X k }. 
Our next step is to prove 

(32) lim n [ E[r 2 nl \X nl = x]dF(x) = 0, Be£. 
First note that 

(33) P{v nl = m\X 1 = x} = P{S n - 1 + l = m}, 

where S n -i is a binomial random variable with parameters n— 1 and F(x). 
Then 



(34) 



n^[r2 1 |X nl =x]=E 



n[l — a(x) 



S„-i+ll2 



(Sn-l + 1) 

To determine the limiting behavior of this expression, we will use the identity 

1 



(35) 



E[f(S n )} 



p(n + 1 



■E[S n+1 f(S n+1 - l)l(5 n+ i > 1)], 



where / : {0, . . . , n} — > 5ft, and is a binomial random variable with param- 
eters n and p. This follows by writing out the expectations. 
Applying (35) to (34) yields 



nE[r 2 a \X n i =x]=E 



[l-a(x) s "] 2 l(S n >l) 



S n F(x) 



The last expression in brackets is bounded by l(S n > l)/(S n F(b)), where 
b = sup{x : x £ B} and it converges to w.p.l since F(b) > and S n — > oo 
w.p.l. Thus, (22) is true by the bounded convergence theorem. 
Next, consider the function 

1 

r(x) = 



a{x)F{x) 



This is bounded on compact sets in E since l/a(x) is, and so the measure \i 
in (27) is bounded on compact sets. It remains to prove (23) in Corollary 5, 
which, in light of (31) and (33), is 



(36) lim / E 

n ^°°JB 



n[l - a(x) Sn ~ 



r(x) 



dF(x)=0, Be£. 



a(x)(S n -i + 1) 

Using (35), we see that the expectation in this integral equals 

■77f„\Sn Sn 



(37) 



r(x)E 



1 -a(xY 



nF(x) 



Hs n > i) 
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The expression inside the expectation is bounded by 1 + 1/F(b), where 
b = sup{x : x £ B} and it converges to because of the classical strong law 
of large numbers n~ 1 S n — > F{x) w.p.l. Then the expectation in (37) con- 
verges to by the bounded convergence theorem. Furthermore, since r(x) 
is bounded on B, the dominated convergence theorem justifies (36). 

This completes the verification of conditions in Corollary 5, which yields 
£ n -i- N as n — ► oo, and hence, Nt — > N as t — > co. □ 

4. Sojourn and departure times of particles. This section describes so- 
journ times of particles in the system defined above. The sojourn times have 
tractable distributions when the arrival times form a renewal or Poisson pro- 
cess. We will also discuss some features of the departure times of particles. 

Consider the particle process {N t : t > 0} described in the preceding sec- 
tion, where the initial state No need not be 0, and Nt may be stationary or 
nonstationary (in the time parameter). We will consider the sojourn time 
W{x) of a "typical" x-particle that arrives at a time t. This sojourn time 
is "conditioned" on having an arrival at a fixed time t with attribute x. 
Conditioning on such events (which may have probability 0) requires the 
use of a family of Palm probabilities Pt (t > 0) for nonstationary systems as 
described in [3]. When a system is stationary, these Palm probabilities are 
independent of t and equal to a single probability measure. Since a formal 
definition of these probabilities is lengthy, we will not include it; the fol- 
lowing discussion can be understood without a working knowledge of Palm 
probabilities. 

For the next result, let Pt denote the Palm probability of the entire system 
conditioned that a particle arrives at time t. In addition to W(x), we will 
consider the spatial sojourn time W(B) of a typical particle whose attribute 
is in the set B £ £. 

Proposition 7 (Sojourn times). Suppose the arrival process A(t) for 
the particle system is a renewal process. Then the Palm distribution of the 
sojourn time W{B) for a particle entering B at time t is 

(38) P t {W(B)>w}= f E[{l-a(x)F(x)) A(w) ]dF(x), w>0. 

JB 

In case A(t) is a Poisson process with rate A, 

(39) P t {W{B)>w}= [ e- Xa ^ x UF(x), 

JB 

and hence, the sojourn time W(x) for an x-particle entering at time t has 
an exponential distribution with rate \a(x)F(x). 
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Proof. By the definitions of the random variables involved, 

(40) Pt{W{x) >w} = P{A(w) <u{x)}, 

where v{x) is the number of new arrivals until the x-particle is deleted. 
Under the rule for particle deletions, an x-particle is deleted by an arbitrary 
arrival with probability a(x)F(x), independently of everything else. Then 
v{x) has a geometric distribution with P\y(x) > n} = (1 — a(x)F(x)) n , and 
u(x) is independent of A(w). Applying these properties to (40), we have 

(41) P t {W(x) >w} = E[(l - a(x)F(x)) A{w) }. 

Then conditioning on x being in the set B yields (38). 

In case A{t) is a Poisson process with rate A, we know that E^s^*)] = 
e -A(i-s) 5 and go ( 4 -g rec j uces to (39)_ □ 

In the preceding result, the sojourn-time distribution of a particle is in- 
dependent of the time at which the particle enters the system. This inde- 
pendence of time is because the sojourn of a particle does not depend on 
any of the system data from the past. It is a rare instance in which a certain 
probability for a nonstationary system is the same as the analogous proba- 
bility for the stationary version of the system. Typically, Palm probabilities 
for a nonstationary system are not equal to those for its stationary version. 
For example, in most input-output systems or queueing networks, customer 
sojourn times depend highly upon quantities such as how many particles (or 
customers) are in the system upon arrival, as well as the total waiting times 
of the particles (customers) present. 

Remark 8. Even though sojourn times of particles in the setting of 
Proposition 7 are identically distributed, they are not independent. This 
follows because the simultaneous deletion of particles triggered by arrivals 
introduce dependencies among the sojourns. 

Proposition 9. The Poisson convergence Nt—> N in Theorem 6 is also 
true when the initial population process Nq is any nonzero point process on 
E. 

Proof. Because the particles do not move about in E, and the distri- 
bution of a point process is determined by its values on compact sets, the 
assertion will follow upon showing that, for any compact set C 6 £, the max- 
imum sojourn time r(C) of the Nq(C) particles in C at time is a finite 
random variable. 

Arguing as in the proof of Proposition 7 and conditioning on Nq{-) and 
A(w), we have 

P{t(C) <w}> E[(l - (1 - b) A{w) ) N ° {C \ 
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where are the locations of the Nq(C) particles and b = inf x6 c a(x)F(x). 
Because A"o is a point process, Nq(C) is finite. Also, A(w) — ► oo w.p.l as 
w — > oo by assumption. Thus, 

P{t(C) < 00} > lim E[(l - (1 - b) A ^) No{C) ] = 1 

w — > OO 

by dominated convergence. This completes the proof. □ 

Another aspect of interest for the particle system is its departure process. 
The departure processes for Jackson and related queueing networks in equi- 
librium are space-time Poisson processes (e.g., see [6]). The departure times 
in the particle system are a little more complicated. The particles depart in 
batches at the arrival times of particles, and each batch is represented by a 
point process on E describing the locations where the particles depart. 

Specifically, the departures are described by the space-time point process 
M on 3? + x E of the form 

M(IxB) = J2HT n eI)M n (B), 

n 

where M n (-) is the point process of batch departures at time T n . Under the 
assumptions for particle deletions, M is a space-time J^-Markov process with 
respect to Tt = a(N s (-) :s <t) [which contains a(A(s) : s < t)]. In particular, 
the increments of M occur at the Poisson arrival times with rate A, and the 
increments M\, M2, ... at these times form a .T^-Markov chain, and M n con- 
ditioned on ^t„_i is a point process with intensity a(x)l(X n > x)Nt ti _ 1 (dx). 

Even when the system is in equilibrium, the increments M n of the depar- 
ture process M are dependent. However, a "typical" batch departure process 
M n is a Poisson process as described below. 

Proposition 10. Suppose the arrival process for the particle system is 
a Poisson process and the system is in equilibrium, and let Mq(-) denote the 
batch departure process conditioned on an arrival to the system at time (as 
defined by the Palm probability). Then the spatial batch- departure process Mq 
under the Palm probability is a Poisson process on E with intensity measure 
a(x)F(x) dfi(x). 

Proof. The Poisson limiting distribution for Nt described by Theo- 
rem 6 is its stationary distribution, since this limit does not depend on the 
distribution of Ao by Proposition 9. So the state of the process in equilibrium 
is given by a Poisson process Af with mean measure [i as in Theorem 6. Then 
under the Palm probability of the process conditioned on an arrival at time 
0, the spatial batch-departure process Mq is a conditional p(x)-thinning of 
N, where p(x) = 1 — a(x)F(x), which is the probability that an x-particle is 
retained in the system. Consequently, Mq is a Poisson process with intensity 
measure a(x)F(x) dfx(x). □ 
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5. Stationary distribution of Markovian particle systems. In the particle 
system we are studying, the probability that an x-particle is deleted by 
an arriving y-particle has the very special form p(x,y) = l(x < y)a(x). A 
natural question is "Are there more general deletion probabilities p(x,y) 
under which Nt has a Poisson limit?" Also, the attribute space E for the 
particles is an interval in 3ft that is linearly ordered. How important is this 
linear order for Poisson limits? 

We will give some insight into these issues with the following results. 
Consider a particle system somewhat like the one we have studied that 
evolves as follows. Particles arrive to the system according to a Poisson 
process A{t) with rate A, and each particle has an attribute in a Polish space 
E. The attributes are i.i.d. with distribution F. Whenever a new particle 
arrives with an attribute y, each particle in the system is subject to being 
deleted; p(x, y) is the probability that an x particle is deleted by the arriving 
y-particle, independently of everything else. The only assumptions we make 
on the probabilities p(x,y) is that J E p(x,y) dF{y) > F-a.e. and that the 
measure 



is finite on compact sets. 

Let Nt(B) denote the number of particles with attribute in B at time t. 
As above, because of the nature of the deletions and the Poisson arrivals, 
{Nt : t > 0} is a Markov jump process. Here is a characterization of the mean 
measure of a point process whose distribution is the stationary distribution 



Proposition 11. Suppose the stationary distribution of the Markov 
process {N t : t > 0} exists and it is the distribution of a point process N 
on E whose mean measure fi is finite on compact sets. Then p, is given by 



Proof. Suppose the Markov process Nt is stationary, which is equiva- 
lent to N = N. Since Nt is a Markov chain subordinated to a Poisson pro- 
cess, its stationary distribution is the same as that of its embedded Markov 
chain N n = N Tn . Therefore, Ni = N. 

Now, Ni = £ + 5x 1 , where X\ is the attribute of the particle that ar- 
rives at time T\, and £ is the point process consisting of those points of 
Nq that survive the deletions at time T\ triggered by the arriving particle. 
Consequently, 



(42) 




of N t . 



(42). 



(43) E[N(B)]=E[N 1 (B)} = E[aB)]+P{X 1 eB}, Be£. 
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Under the assumptions on particle deletions, we know that when X\ = x%, 
the process £ is a p(x, xi)-thinning of iVo = N. Then conditioning on X±, 



E[£(B)] = E / [l-p(x,X 1 )]dN(x 

= H(B)- I E\p{x,X x )]dii{x). 

JB 

Substituting this in (43) yields 

/ E\p(x,X 1 )]dn(x) = P{X 1 eB}. 
Jb 

Thus, fi is given by (42). □ 



We will now consider conditions under which the stationary distribution 
of the Markov process Nt is that of a (spatial) Poisson process on E. We 
already have one example. 

Remark 12. Suppose Nt is the particle process as in Section 3 with a 
Poisson arrival process with rate A. The stationary distribution of Nt is that 
of a Poisson process iV with mean measure given by (27). This follows since 

Nt — ► N by Theorem 6, and this limit is independent of the distribution of 
Nq by Proposition 9. 



Proposition 13. The stationary distribution of the Markov process 
{Nt : t > 0} is that of a Poisson process on E with mean measure (42) if 
and only if 

(44) / ei^^W^-ZW dF( Xl ) = 1, / € C+ 
Je 

where h(x,xi) =p(x,x 1 )/ J E p(x,y)dF(y). 

Proof. Suppose A^o is a Poisson process on E with mean measure \i 
given by (42). Then the stationary distribution of the Markov process Nt is 

that of the Poisson process iVo if and only if N\ = N$. The latter, in terms 
of Laplace functionals, is equivalent to 

(45) E[e-^f] = E[e- N °f], f £ C+. 

Then to prove the proposition, it suffices to show (45) is equivalent to 
(44). Since iVo is a Poisson process, 

E[e- N °f] = eS* ( - 1 - e ~ I{ ' ) M x \ feC+. 
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Using the representation Ni = £ + 5x ± , from the preceding proof, and con- 
ditioning on Xi, we have 

(46) E[e-* lf ] = I E[e-tf\X 1 = x l ]e- f[xi) F{dx l ). 

By the assumptions on particle deletions, we know that when X\ = Xi, 
the process £ is a p(x, rri)-thinning of the Poisson process Nq. Then the well- 
known theorem on thinning of Poisson processes implies that £, conditioned 
on X\ = x\, is a Poisson process with mean measure [1 — p{x,x\)\ dfi(x). 
Therefore, 

]QgE[e-* f \X 1 = xi] = - I (l-e- /(x) )[l-p(x,X!)]^(x) 

JE 

= logE[e- N °f]+g(x 1 ), 

where 

g(xi)= / {l-e- f{x) )p(x,x 1 )dn(x). 
Je 

Substituting this expression in (46), we obtain 

(47) E[e-* lf ]=E[e- Nof ] f e~ [fi - Xl) - 9{xi)] F(d Xl ). 

Je 

Therefore, (45) is equivalent to 

Je 

But this expression, with dfj,(x) = dF(x)/ J E p(x,y) dF(y) from (42), is equiv- 
alent to (44). □ 

Example 14. The particle process in Section 3 satisfies condition (44). 
Hence, the stationary distribution of Nt is that of a Poisson process, which 
we had already confirmed in Remark 12. This observation follows, since in 
this case, 

p(x,y) = l(x < y)o(x), h(x,xi) = l(x < x\)/F{x), 
and (44) reduces to 

(48) [ r^W^!^) = 1. 

Here one uses "w~r = ~~ logi ? (xi). Then clearly (48), with a change of 
variable, equals J °° e~ u du = 1, which is true. 
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Example 15. Counterexample. Consider a particle system that evolves 
as the process described in Section 2, except that an arrival to y will delete 
an x-particle with probability p(x,y) = a(x), independent of y. In this case 
h(x, x\) = 1 and (44) reduces to / E e~^ x > dF(x) = 1, which is not true for all 
/ £ Cft- Hence, the stationary distribution of Nt cannot be that of a Poisson 
process. 
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